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Abstract 

We propose a simple quantum mechanical model describing the time de- 
pendent diffusion current between two fermion reservoirs that were initially 
disconnected and characterized by different densities or chemical potentials. 
The exact, analytical solution of the model yields the transient behavior of 
the coupled fermion systems evolving to a final steady state, whereas the 
long-time behavior is determined by a power law rather than by exponen- 
tial decay. Similar results are obtained for the entropy production which is 
proportional to the diffusion current. 

1 Introduction 

Key ingredients to the investigation of non-equilibrium features of physical sys- 
tems, time irreversibility, energy dissipation and entropy production [1, 2, 3, 4, 5] 
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can be treated analytically only in a few cases. Providing a well-known example 
of the latter, the Caldeira-Leggett model [6] along with all its extensions [7, 8, 9] 
essentially leans on the quadratic interaction imposed to couple a huge number 
of harmonic oscillators. As such, the resulting equations of motion governing the 
dynamics of any relevant observables remain classical, in spite of the quantum me- 
chanical nature of the density matrix, potentially representing the initial equilib- 
rium conditions. In this light, it is tempting to conceive an exactly solvable model 
that reveals a signature of quantum dynamics while exhibiting at least one of the 
typical non-equilibrium features such as irreversibility. The model presented be- 
low describes quantum diffusion as the basic mechanism underlying the particle 
exchange between two fermion systems that are brought into contact at some time 
instant t = 0, while having been separated at all previous times. Section 2 intro- 
duces the model in terms of its Hamiltonian and discusses the time evolution in 
the framework of the Heisenberg picture. In section 3 the Heisenberg equations 
of motion for the fermion creation and annihilation operators are explicitly solved 
and the time dependent obervables, including the diffusion current, are extracted 
by averaging the corresponding Heisenberg operators with the initial density ma- 
trix characterizing the initial equilibrium of the uncoupled fermion reservoirs. As 
a result, the time dependent fermion densities are calculated explicitly and the re- 
sults are discussed in terms of the occurrence of time irreversibility. In particular, 
a direct relation between the diffusion current and the entropy production is ex- 
ploited to study the time dependence of the latter. Finally, a conclusion is drawn 
in section 4. 



2 Two coupled fermion systems: Heisenberg equa- 
tions of motion 

We consider two non-interacting fermion gases establishing two-dimensional reser- 
voirs in thermal equilibrium that are completely disconnected at all times preced- 
ing t = 0, while being confined to two rectangular areas with sizes L x and L y . The 
uncoupled reservoirs are at the same temperature, however with different initial 
densities n w and n 2 o- 

The Hamiltonian of the disconnected 2DEGs can be written as 

H = H i + H 02 = J] 6 k (cj k c lk + <4 k c 2k ) , (1) 

k 
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where c lk ,c\ k and c 2 k,C2 k respectively denote the fermion creation and annihila- 
tion operators for both subsystems and e k are the single-particle eigenergies. For 
the sake of simplicity, spin degrees of freedom have been ignored whereas the 
eigenenergies are taken to be parabolic, i.e. e k = fi 2 k 2 /2m with m denoting the 
fermion effective mass. As of t = 0, the two reservoirs get connected by turning 
on the perturbation 

A ' = Z Ak ( c ^ c 2k + 4 c ik)- ( 2 ) 

k 

which initiates the exchange of fermions. A k is the transition rate specifying the 
number of fermions with a given wave vector k being transferred from one reser- 
voir to the other. Throughout this paper A k is considered real and even in k. The 
total Hamiltonian now reads 

H = H O + 0(OH'. (3) 

The time dependence of the observables characterizing the diffusion process will 
be extracted from the time dependent Heisenberg operators ci k (t), c\ k (t), C2k(t), 
cl k (t). The latter are known to obey the Heisenberg equations of motion 



dc jk (t) 

iti — — = [c jk (t),H], ift 
dt 

which, in the present case, reduce to 



d ^kW 
dt 



= [c] k (t)M j = 1,2. 



ift dC ' k ^ = cikO) + A k c 2 k(0. 
dt 

dC2k(0 

\fi —7— = A k cikO) + e k c 2k (0, 
dt 



(4) 

(5) 
(6) 



because the perturbation (Eq. 2) does not mix momenta during a transfer process. 
Consequently, we may analytically solve the equations of motion and explicitly 
present the solutions in terms of the Heisenberg operators at t = 0, thus yielding 
separately the time evolution in the form of 2x2 unitary matrix for each k, 









= exp| 


,C 2 k(0, 
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-i sin 
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cos 



H 



(c lk (0) 

C2k(0)j 



(7) 
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In particular, the operators the average of which provides the occupation numbers 
of the first subsystem are given by 



cL(t) c lk (0 = cos 2 



(A*fi 



cj.(0)c lk (0) + sin 2 



(A k t 



4(0)c 2k (0) 



(A k t 



- 1COS 



(A k t 



sin 



[4(0)c 2k (0)-4(0)c lk (0)] (8) 



whereas a similar expression follows for c 2k (t) c 2k (t). As all time dependence is 
in the operators, the Heisenberg picture only requires the explicit knowledge of 
the initial density matrix p () to calculate any expectation values at arbitrary later 
times. In general, a time dependent observable Y, representing the time dependent 
ensemble average of a Heisenberg operator Y(t) is to be extracted according to the 
recipe [10] 

Y t = (Y(t)) = Tr[p Y(t)]. (9) 

As both fermion reservoirs are assumed to be uncoupled at t = and before, 
the initial density matrix p reduces to the product of the two reservoir density 
matrices 

p = ^-exp(-jS(Hoi - mNi)) exp(-^(H 02 - u. 2 N 2 )) 

= y HI ex P ("^ ^ 6k ~ ^ c ik c ik) exp (-p (e k - \i 2 ) 4 c 2k) . (10) 

where ft = l/k B T, Z is the partition function and Ui,U2 denote the reservoir 
chemical potentials compatible with the initial densities n w and n 20 . The number 
operators counting the fermions in the distinghuised reservoirs are given by 



c j k c jk, 



7 = 1,2. 



(ID 



Obviously, the initial density matrix (10) translates to the Fermi-Dirac distribution 
ruling the initial occupation numbers: 

1 



(cl k (0)c jk (0)) = F j (e k )6 l j, Fj(E) = 



1+exp (j6(£-^))' 



7,/= 1,2. (12) 



Combining Eqs. (8) and (12), we arrive at 



(A k t 




'A k f 


[nj 


F x (6k) + sin 2 
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F 2 (e k ), 



(13) 



In particular, we obtain the generic expression for the time dependent fermion 
density in the first reservoir from 



m(t) 



1 T (mt)) = ^j-Yj(c\ k (t)c lk (t)) 

x y k 



L X Ly 



1.1,. Zj 



Lj(Ly 

which can be rewritten as 





( M 7 


fA k ^ 




cos 2 


— 1 + sin 2 


[ *J 



(14) 



1 

2 



= - (»io+«2o)+ ' , y 



'2A k t 



COS 



ft 



[Fi(^)-F 2 (^)] (15) 



3 Solution and discussion 

Defining the diffusion current density J D (t) to be the rate of change of the fermion 
density in the first reservoir, we may obtain the corresponding expression for the 
present model by simply taking the time derivative of Eq. (15), 



J D (t) = 



d«i(0 
dt 



1 



— 2j A k sin 



( 2A k ^ 



[F 1 (e k )-F 2 (e k )] 



(16) 



It is worth noting that the entropy production is unambigously determined by the 
total diffusion current 7 D = J D L x L y = d(Ni(0) / dt, the initial difference of the 
chemical potentials \i\-\i2 and the common, ambient temperature T. Starting from 
the time dependent density matrix p ; satisfying the quantum Liouville equation 



»£-[<U], 

we first consider the relative entropy [1, 2, 3, 4, 5] 

S [p t , Pol = -^Tr [p, (In p - In p f )] . 



(17) 



(18) 



From Eq. (17) one may easily derive the well-known result Tr(p ; In p f ) = Tr(p In p ) 
which is the mathematical statement of the Shannon entropy failing to evolve in 
time. Substitution into Eq. (18) yields 



S [pt, Po] = ~k B Tr [(pt - po) In p ] , 



(19) 
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whereas the factorized, initial density matrix in Eq. (10) leads to 

Tr (p t In p ) = -ySTr [p, (h - ^Nj - |iiN 2 )] . (20) 

Next, we use p, = exp (-iHt/fr) p exp (iHt/fr) which provides the formal solution 
to the Liouville equation, to eliminate the time dependence of the density matrix 
and recast Eq. (20) in the Heisenberg picture: 

Tr(p t lnp ) = -/?Tr[po(H (0 - pjN^) - u. x N 2 (0)] . (21) 

Using Eq. (21) and bearing in mind that H is a constant of motion because the 
diffusion process is not accompanied by energy exchange, we obtain 

S [p t , Pol = ~ J] & (<N/0> - (Ny(0)))] . (22) 

7=1,2 

Eq. (22) reveals that any entropy change is to be considered configurational in 
the sense that it is exclusively due to the redistribution of all fermions over both 
reservoir areas. Emerging as the time derivative of the relative entropy, the entropy 
production can be directly derived from Eq. (22): 



dS(t) 



=--y 

dt T k 



d<N,<0> 



(23) 



Finally, taking into account that the total number of fermions is a constant of 
motion as well, we arrive at 

^ = -I(m-^ 2 )/ D (0. (24) 

Hence, the entropy production is instantly proportional to the diffusion current 
while the sign of the latter depends on which of either reservoirs is being emp- 
tied into the other. In order to obtain more detailed results or numerical output 
based on Eq. (15), it is paramount to specify the functional dependence of A k on 
the wave vector k as well as to identify the range of allowed wave vectors. A 
simple and straighforward model favoring the transfer of high-energy fermions 
corresponds to A k being chosen proportional to e k , say 

A k = A k = l €k , l>0 (25) 
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and will be further explored as a benchmark. Moreover, it follows from Eq. (15) 
emerging as a sum of trigonometric oscillations that the time evolution of the 
present diffusion model is bound to be reversible as long as the allowed wave 
vectors run through a finite, countable set. Phrased otherwise, time irreversibility 
is expected to occur only if the thermodynamic limit is taken with L x ,L y , (Ni), 
(N 2 ) —> oo, the densities n\ and n 2 however remaining finite. Traditionally, this 
encompasses the conversion of the summation over k into a 2D integration ac- 
cording to 

Carrying out the latter and transforming to the "polar" coordinates (Ve, (p) ac- 
cording to 



V2m£ , V2m£ 

— — cos (p, ky = — — 

n ' n 



k x = — - — cos <p, k y = — - — sin 0, E > 0, < <f> < In, (27) 



we obtain 



it) = \ (mo + «2o) + ^ Jd 2 k cos [Fxte) - F 2 (e k )] 



1 m f°° C 2n (2XEt. 

= ^ ("io + «2o) + J dE ^ dcf> cos I — ) - F 2 (£)] 



1 m f°° /2XEA 

= ^ ("io + «2o) + ^5 J d£ cos I — J [Fx(£) - F 2 (£)] . (28) 

Clearly, the damped oscillations are driven by the initial chemical potential dif- 
ference \ii - \i 2 - Note that the 2D geometry of the two reservoirs enables us to 
extract analytically the relation between \i 2 and the corresponding initial den- 
sities from 

exp , ' -1 



[ij = k B T\n 



J =1,2. (29) 



mk B T 

Any further analytical treatment of Eq. (28) amounts to either a series expansion of 
the Fermi-Dirac distribution function when it comes to keep track of the temper- 
ature dependence, or to studying the extreme quantum limit at zero temperature. 
In the latter case, the Fermi-Dirac functions reduce to step functions respectively 
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imposing integration boundaries at E = \ii > and E = \x 2 > 0, yielding 
1 m 1 [ . /2Xii 1 A . (2X\i 2 t 

= X («10 + «2o) + TT^-T • - Sln - Sln 



2 f I ft / I ft 



(30) 



On the other hand, in the classical limit corresponding to k B T » nft 2 nj Q /m and 
(J,i,(i2 < we may replace the Fermi-Dirac functions with bare exponentials to 
arrive at 

m(f) = 5 («io + n 20 ) + £^ (e^> - e^ l2 )j^d£ cos e** (31) 

^ ( „ 10+ „ 20) + ^ (e ^_ e ^)__i__ 

= « ("10 + «20) + X ("10 - «2o) 



2 V 1U zu/ 2 V 1U /3 2 ti 2 + 4A 2 t 2 

This time, the evolution of n\{t) towards its steady-state value n s = 1 /2(n\o + n 2 o) 
is governed by a r 2 power law, while no oscillations are observed during the tran- 
sient regime. Note that, whenever irreversibility is established, the evolution from 
the original equilibrium state to the final steady state corresponds to an entropy 
increase because n l0 , the original density in reservoir 1, decreases (increases) to- 
wards n s provided that > (<) [i 2 , as can also be seen from the entropy change 
per particle, 

* ,. S[p t , po] U.1-U.2 n w -n 2Q 

As = lim — = — . (32) 

t^oo (Ni(0) + <N 2 (0> T + "20 

4 Conclusion 

The quantum diffusion model presented in this paper can be solved analytically 
to obtain a closed-form integral representation of the time dependent observables, 
such as the fermion density in the reservoirs. The long-time behavior is governed 
by a power law t~ a where a = 1 and a = 2 respectively characterize the extreme 
quantum regime and the classical regime. The time dependent entropy production 
is proportional to the diffusion current while the irreversibility of the latter in the 
case of the thermodynamic limit corresponds to a net increase of the (configura- 
tional) entropy when the steady- state is attained. 
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